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Abstract- As  a  powerful  tool,  the  symmetric  phase-only  matched 
filter  (SPOMF)  has  been  shown  to  yield  superior  performance 
over  the  conventional  correlator  and  is  widely  used  in  image 
registration  and  recognition.  In  this  paper,  we  investigate  the 
use  of  this  SPOMF  for  processing  GNSS  signals.  This  extension 
is  compatible  with  our  frequency-domain  software  GNSS 
receiver  architecture  in  which  both  the  incoming  signal  and 
replica  spectra  are  available  for  the  SPOMF  implementation 
versus  the  conventional  correlator. 

The  use  of  phase-only  information  is  equivalent  to  equalizing 
the  magnitude  spectrum  in  contrast  to  the  original  spectrum 
that  tapers  off  according  to  a  s/nc-function.  This  tends  to 
accentuate  the  high  frequency  components  corresponding  to 
edges  or  transitions  in  the  signals.  As  such,  the  SPOMF 
produces  a  much  sharper  peak  (ideally  a  Dirac  delta  function) 
that  is  more  accurate  in  timing  and  less  sensitive  to  multipath. 
In  addition,  the  same  operation  is  applicable  to  both  a  binary 
phase  shift  keying  (BPSK)  modulated  signal  such  as  the  GPS 
C/A-code  and  P-code  and  a  binary  offset  carrier  (BOC) 
modulation  such  as  the  GPS  M-code  and  Galileo  codes.  More 
importantly,  it  only  has  a  single  matching  peak  regardless  of 
which  modulation  code  is  being  used. 

In  this  paper,  the  SPOMF  is  introduced  within  the 
framework  of  a  generalized  frequency-domain  correlator 
(GFDC)  for  GNSS  signals.  The  salient  features  of  SPOMF  as 
well  as  its  application  to  BPSK  and  BOC  signals  are  illustrated 
with  simulation  examples. 

I.  Introduction 

Correlation  is  a  critical  operation  in  GPS  receivers.  As  a 
spread-spectrum  signal,  the  received  GPS  signal  is  below  the 
thermal  noise.  Through  despreading  integration,  correlation 
provides  the  processing  gain  necessary  to  detect  a  signal  from 
noise  and  by  doing  so,  it  identifies  which  satellite  the 
acquired  signal  is  coming  from.  Correlations  with  code  and 
carrier  replicas  at  different  code  phases  and  carrier  phases 
provide  code  delay  errors  and  carrier  phase/frequency  errors 
that  are  used  as  input  to  code  and  carrier  tracking  loops, 
respectively.  The  accumulated  correlations  are  further  used 
for  navigation  data  bit  sync  and  demodulation.  These  and 
other  pieces  of  information  lead  to  GPS  observables  and 
satellite  ephemeredes  for  the  ultimate  timing  and  position 
fixing  [5,  7,  10,  12]. 

For  a  maximum-length  pseudo-random  number  (PRN) 
code,  its  correlation  function  is  ideally  an  equilateral  triangle 
with  its  base  width  being  ±TC  where  Tc  =  l/fc  is  the  chip 


duration  and  fc  is  the  chipping  rate.  In  the  acquisition  mode, 
the  search  step  in  code  phase  is  typically  chosen  as  At  =  T</2 
and  this  produces  the  worst  signal  loss  of  2.5  dB.  It  is 
therefore  desirable  to  have  a  correlation  function  with  a  wide 
base  in  the  acquisition  mode  so  that  a  minimal  number  of 
steps  are  needed  to  cover  a  given  interval  of  time  uncertainty. 
That  was  one  reason  for  which  the  GPS  C/A-code  (i.e.,  the 
coarse  acquisition  code)  was  originally  designed  so  as  to 
assist  the  precision  P-code. 

In  the  tracking  mode,  however,  it  is  desired  to  have  a 
correlation  function  with  its  base  as  narrow  as  possible.  Such 
a  sharp  correlation  function  not  only  provides  more  precise 
timing  but  also  results  in  a  delay  estimate  that  is  less  sensitive 
to  multipath.  This  is  due  to  the  fact  that  only  very  closely 
spaced  multipath  signals  (<  1.5TC )  can  contribute  to  the  delay 
estimation  errors  and  at  much  closer  ranges  the  errors  are 
insignificant  as  compared  to  thermal  noise. 

Clearly  different  requirements  are  imposed  on  correlation 
function  for  acquisition  and  tracking,  which  are  in  conflict.  A 
code  that  may  satisfy  both  requirements  at  the  same  time  is 
the  binary  offset  carrier  (BOC)  modulation  as  used  by  the 
new  M-code  and  LIC-code  (under  development)  as  well  as 
the  European  Galileo  codes.  By  placing  the  signal  power 
toward  the  edges  of  the  frequency  band,  it  was  designed  to 
co-exist  with  binary  phase  shift  keying  (BPSK)-codes  already 
in  the  same  band  with  reduced  interference  since  a  BPSK 
code  power  is  concentrated  around  the  band  center.  The  PRN 
code  of  a  BOC  modulation  has  a  correlation  function  with  a 
large  base  that  can  be  obtained  with  a  single  sideband  (either 
upper  or  lower)  [16,  17].  Together  with  the  offsetting  square 
wave,  the  composite  correlation  of  a  BOC  code  has  a  refined 
mainlobe.  However,  it  also  has  numerous  nulls  and  sidelobes. 
Although  the  mainlobe  is  narrow,  the  sidelobes  are  not 
substantially  smaller.  Without  special  hardware  and  software, 
a  receiver  runs  the  risk  of  being  trapped  in  nulls  (i.e.,  missing 
detection)  or  locking  onto  a  sidelobe  (i.e.,  biased 
measurements)  [1,  3]. 

In  this  paper,  we  set  forth  a  generalized  frequency-domain 
correlator  (GFDC)  to  satisfy  the  correlation  function 
requirements  in  both  the  acquisition  and  tracking  modes  for 
BPSK-  and  BOC-types  of  codes  alike.  The  refined  correlation 
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function  is  obtained  by  implementing  the  symmetric  phase- 
only  matched  filter  (SPOMF)  in  the  tracking  mode  while  a 
conventional  correlator  is  implemented  in  the  acquisition 
mode  (with  single  sideband  segmentation  for  the  BOC  codes 
for  instance  [16]). 

The  SPOMF  is  widely  used  in  image  registration  and 
recognition  [2,  9,  13]  and  its  extension  to  GNSS  signal 
processing  is  compatible  with  our  frequency-domain  GNSS 
software  receiver  architecture  [15]  in  which  both  the 
incoming  signal  and  replica  spectra  are  available  for  flexible 
implementation.  In  this  paper,  the  SPOMF  is  introduced 
within  the  framework  of  GFDC  for  GNSS  signals.  The  salient 
features  of  SPOMF  as  well  as  its  application  to  BPSK  and 
BOC  signals  are  illustrated  with  simulation  examples. 

n.  Generalized  Frequency-Domain  Correlator  (GFDC) 

We  will  introduce  the  symmetric  phase-only  matched  filter 
(SPOMF)  as  a  particular  implementation  of  the  generalized 
frequency-domain  correlator  (GFDC),  as  part  of  a  frequency- 
domain  software  GNSS  receiver.  However,  it  can  be  used 
standalone  or  as  a  replacement  of  conventional  correlators  in 
conventional  GPS  receivers. 

The  inherent  flexibility  of  a  software  GPS  receiver  allows 
it  to  adopt  a  frequency-domain  baseband  processor  for  which 
a  GFDC  is  shown  in  Fig.  1.  The  use  of  a  GFDC  has  several 
advantages.  First,  the  same  operation  is  applicable,  without 
any  other  changes  except  for  the  replica  code,  to  both  the 
BPSK  codes  such  as  GPS  C/A-code  and  the  BOC 
modulations  such  as  M-code.  Second,  it  can  switch 
seamlessly  among  different  types  of  correlation/matching 
operations  with  the  wide  base  for  acquisition  and  narrow  base 
for  tracking.  Third,  the  generalized  correlation  has  a  sharp 
peak  when  the  SPOMF  is  in  use.  The  reduced  base  width 
makes  it  more  accurate  in  timing  and  less  sensitive  to 
multipath. 

This  GFDC  differs  from  the  generalized  cross  correlator 
(GCC)  that  has  been  used  for  radar  and  sonar  signal 
processing  for  delay  estimation  [4,  6].  The  GCC  is  closely 
related  to  the  coherence,  a  complex  quantity  that  is  the  cross¬ 
power  spectral  density  between  two  random  processes 
divided  by  the  product  of  their  auto  power  spectral  densities. 
The  design  goal  of  the  GCC  is  to  produce  the  highest  signal- 
to-noise  ratio  (SNR)  correlation  peaks  in  the  presence  of 
noise.  The  optimum  filter  is  designed  with  the  spectral 
characteristics  of  the  noise  assumed  to  be  known.  However, 
some  of  these  filters  alter  the  phase  of  the  input  functions, 
which  could  bias  the  final  output  [13]. 

As  shown  in  Fig.  1,  the  GFDC  accepts  the  incoming  signal 
s(t)  and  the  code  replica  r(t)  and  produces  the  correlation 
function  c(t)  between  the  two.  This  is  the  same  input-output 
behavior  as  other  correlators.  However,  what  makes  it 
different  from  others  is  the  spectrum  filtering  U(f),  Iff),  and 
W(f)  applied  along  the  signal  processing  chain  and  the 
feedback  paths  inserted  in  the  frequency  domain. 

Clearly,  when  U(f)  =  1,  V(f)  =  1,  and  W(f)  =  1  (i.e., 
without  any  spectrum  manipulation  and  feedback),  it 
becomes  a  simple  straight  FFT-implemented  correlation. 


Various  spectrum  filtering  can  be  applied  that  makes  the 
GFDC  design  versatile,  which  is  discussed  below. 


Fig.  1 .  Generalized  Frequency-Domain  Correlator  (GFDC)  Architecture 

A.  Incoming  Signal  Spectrum  Filtering  and  Feedback 

A  number  of  spectrum  filtering  techniques  applicable  to 
the  incoming  signal  are  described  below.  They  do  not  serve 
the  same  purposes  and  not  all  are  applicable  at  the  same  time. 

(a)  Spectrum  Excision  of  Narrowband  Interference.  Being 

spread-spectrum,  the  GPS  signal  is  below  the  thermal  noise. 
Any  spikes  in  the  signal  spectrum  are  attributed  to 
interference  and  as  such,  the  spectrum  values  at  those 

frequency  bins  can  be  removed  and  replaced  with  zeros.  This 
requires  continuous  monitoring  of  the  spectrum  and  power- 
detection  against  a  pre-defined  threshold  [19].  The  operation 
of  excision  can  be  viewed  as  passing  the  signal  through  a 
notch  filter  in  the  time  domain  or  equivalently  as  multiplying 
the  spectrum  with  zero  at  those  frequency  bins  of 

interference.  Assume  there  are  N  frequency  bins.  The  N 

spectrum  complex  values  can  be  put  into  a  vector  as  S  = 
[S(f),  f  =  0,  1,  ...,  N-1]t,  where  the  superscript  T  stands  for 
transpose.  To  remove  an  interference  component  at  the  kth 
frequency  bin,  the  following  operation  is  applied: 

Sk=ZkS  (la) 

Zt=diag([ iM  0  Hj)  (lb) 

where  diag(v)  stands  for  a  diagonal  matrix  with  its  diagonal 
elements  specified  by  y  and  7„  is  a  vector  of  n  ones  (1  ’s).  The 
operation  can  be  repeated  for  all  frequency  bins  of  interest.  It 
is  implied  that  the  same  operation  is  applied  to  both 

corresponding  positive  and  negative  frequency  bins  at  the 
same  time.  Zone-zeroing  and  individual  excision  are  two 
popular  ways  to  apply  this  spectral  filtering. 

(b)  Spectral  Filtering  to  Reduce  Additive  Noise.  To  restore 
the  signals  under  additive  noise,  the  following  spectral  filter 
can  be  applied: 

mf)j  p*(/)  T  (2a) 

where  Pff)  and  P„(f)  represent  the  power  spectrum  of  the 
signal  and  the  additive  noise,  respectively.  When  a  =  (3  =  1, 
Eq.  (2)  corresponds  to  Wiener  filtering.  When  a  =  1  and  f  = 
Z2,  it  applies  power  spectrum  filtering.  It  is  clear  from  Eq.  (2) 
that  the  filter  is  non-causal  with  an  even,  real  frequency 
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response  and  consequently  the  phase  of  the  filter  is  zero.  This 
filter  does  not  affect  the  signal  phase  but  does  modify  the 
magnitude. 

When  the  signal  plus  noise  power  spectrum  is  estimated 
from  the  observed  samples,  denoted  by  p  (y ) ,  Eq.  (2a)  can  be 
modified  as: 

U(f)=  .  1  (2b) 

PJf) 


frequency  uncertainty  interval.  However,  in  the  tracking 
mode,  the  estimated  Doppler  and  its  uncertainty  will  reduce  it 
to  a  small  interval,  say,  from  do-1  to  do+1  where  do  is  the 
closest  bin  to  the  estimated  frequency. 

The  feedback  shown  in  Fig.  1  indicates  the  need  to  repeat 
the  operation  for  each  frequency  bin,  producing  the  delay- 
Doppler  map  of  generalized  complex  correlations. 

(e)  Spectrum  Windowing.  Two  possible  window  functions 
are: 


This  spectral  filter  can  be  used  to  suppress  narrowband 
interference  [11]. 

(c)  Spectrum  Segmentation  of  Multiple  Codes.  In  the  case 
where  each  GPS  signal  band  (LI,  L2,  or  L5)  is  processed,  the 
spectrum  for  individual  codes  need  to  be  extracted  from  the 
frill  band  spectrum  [16].  The  operation  can  be  viewed  as 
passing  the  signal  through  a  bandpass  filter  in  the  time 
domain  or  equivalently  as  selecting  the  spectrum  at  those 
frequency  bins  of  interest.  Both  are  linear  operations.  Assume 
we  want  to  segment  the  spectrum  at  frequency  bins  from  i  to 
j.  The  following  matrix  multiplication  implements  an  ideal 
bandpass  filtering: 

S{ =G/S  (3a) 

Gj=diagmi -+i,Qv_,])  (3b) 

where  0„  is  a  vector  of  n  zeros  (0’s). 

(d)  Spectrum  Translation  for  Residual  Doppler  Removal 
with  Feedback.  The  residual  Doppler  in  the  incoming  signal 
appears  as  a  multiplicative  sine  or  cosine  term  to  the  code 
sequence,  thus  introducing  a  phase  change  from  sample  to 
sample.  It  is  typically  removed  by  multiplying  the  incoming 
signal  with  a  carrier  replica  in  the  form  of  complex 
exponential  at  the  desired  Doppler  frequency.  This  time- 
domain  phase  rotation  is  equivalent  to  spectrum  translation  in 
the  frequency  domain  [19].  To  remove  a  residual  Doppler  of 
±dAf  Hz  where  Af  is  the  frequency  resolution  (i.e.,  the  width 
of  each  frequency  bin),  the  spectrum  is  down  (up)-translated 
by  d  bins. 

Sd=TdS  (4a) 

Td=0  =  I  (4b) 

0  0,_2  i  oN_f 

0  0V_3  0  1  (4c) 

1  0  0V_3  0 

—  d-2  1  0 

where  the  first  non-zero  element  in  the  first  row  is  at  the  dlh 
column  when  d  >  0  and  at  the  ( N+d)th  column  when  d  <  0, 
with  N  being  the  number  of  frequency  bins.  Since  the  matrix 
multiplication  in  Eq.  (4a)  is  equivalent  to  index  permutation, 
the  practical  implementation  resorts  to  circularly  shifted 
indexing  of  the  array. 

Without  knowing  the  residual  Doppler  frequency,  the 
initial  search  will  select  a  large  value  for  d  so  as  to  cover  the 


(5a) 

(5b) 


where  Uo  is  a  real  constant.  Since  the  GPS  signal  spectrum 
falls  off  with  increasing  frequency,  the  inverse  magnitude 
window  frinction  of  Eq.  (5a)  acts  as  a  high-pass  filter,  which 
tends  to  emphasize  the  edge  information  without  affecting  the 
phase  information.  Since  high  frequency  information 
decorrelates  quickly,  it  helps  to  produce  a  very  sharp 
correlation  peak.  This  operation  equalizes  the  magnitude 
spectrum,  thus  keeping  the  phase-only  information. 

B.  Code  Replica  Spectrum  Filtering 

Instead  of  down-converting  the  incoming  signal  from  IF  to 
baseband,  it  is  possible  to  up-convert  the  code  replica  from 
baseband  to  IF  so  as  to  catch  up  with  the  incoming  signal  by 
accounting  for  the  unknown  Doppler  shift.  This  operation  is 
similar  to  the  spectrum  translation  described  above. 

(a)  Spectrum  Windowing.  Window  functions  can  also  be 
applied  to  the  code  replica  spectrum  as  shown  in  Fig.  1.  Two 
possible  window  functions  are: 


V(fY 


K 

Rif) 


(6a) 

(6b) 


where  Vo  is  a  real  constant.  It  has  the  same  effects  as  on  the 
incoming  signal  spectrum  described  above  for  Eq.  (5). 

(b)  Amplitude  Compensated  Filtering.  A  nonlinear  filter 
can  be  used  to  attenuate  any  dc  information  near  the  origin  on 
the  frequency  axis: 


Rq 

I2 

.!«(/)! 


I  «(/)!>  «» 

l*(/)Mo 
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where  Ro  is  a  positive  constant  selected  to  maximize  the 
signal  to  noise  ratio  [8]. 


C.  Correlation  Spectrum  Filtering  and  Feedback 

A  window  frinction  can  be  applied  to  the  correlation 
spectrum.  This  will  shape  the  correlation  frinction  after  the 
inverse  Fourier  transform  is  taken.  Design  goals  include  (1) 
reducing  the  sidelobe  level,  (2)  maximizing  the  correlation 
peak,  and  (3)  narrowing  the  mainlobe.  These  goals  typically 
are  not  compatible  to  each  other.  For  example,  reducing  the 
sidelobe  level  comes  at  a  price  of  enlarging  the  mainlobe. 
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Another  useful  operation  is  to  evaluate  the  correlation 
function  at  a  desired  code  lag.  In  the  acquisition  mode,  since 
the  code  phase  is  unknown,  a  large  interval  of  code  phase 
(timing)  uncertainty  is  searched.  However,  once  in  the 
tracking  mode,  the  estimated  code  phase  and  its  uncertainty 
can  be  used  to  determine  code  lags  where  the  correlation 
function  needs  to  be  evaluated  [14].  The  feedback  shown  in 
Fig.  1  indicates  the  number  and  location  of  the  code  lags  to 
evaluate. 

III.  Symmetric  Phase-Only  Matched  Filter  (SPOMF)  Features 

Given  the  possible  ways  to  filter  the  incoming  signal, 
replica,  and  correlation  spectra  jointly  or  independently,  we 
consider  four  cases  of  particular  interest  to  the  GFDC 
implementation. 

(a)  Correlation  Function.  In  the  first  case  where  U(f)  = 
V(f)  =  W(f)  =  1,  the  generalized  correlation  spectrum 
becomes: 

C(f)  =  S(f)R*(f)  (8a) 

c(t)  =  IFFT{C(f)}  (8b) 

where  the  superscription  *  stands  for  complex  conjugate. 

Eq.  (8a)  is  the  conventional  correlation  spectrum  and  Eq. 
(8b)  is  just  the  FFT-implemented  cross-correlation. 

(b)  Signal  Channel  Transfer  Function/Impulse  Response. 
In  the  second  case  where  U(f)  =  V(f)  =  l/\R(f)\  and  W(f)  =  1, 
the  generalized  correlation  spectrum  becomes: 

C(f)_S(f)R'(f)_S(f)  (9a) 

!«(/)  I2  Rif) 

h(t)  =  I FFT { C(f)}  (9b) 

Eq.  (9a)  is  the  GPS  signal  channel  transfer  function  and 
Eq.  (9b)  is  the  corresponding  impulse  response. 

Although  V(f)  =  l/\R(f)\  effectively  extracts  the  phase-only 
information  from  the  replica  R(f),  U(f)  may  fail  to  do  so 
because  S(f)  generally  differs  from  R(f)  due  to  residual 
Doppler  and  noise.  By  consequence,  some  spectral  spikes 
show  up  near  multiples  of  chipping  rate  where  noise  is  known 
to  dominate.  To  suppress  such  noise  amplification  effects, 
spectral  filtering  such  as  zone-zeroing  and  individual  excision 
can  be  applied  [17,  18]. 

(c)  Phase-Only  Matched  Filter  (POMF).  In  the  third  case 
where  U(f)  =  1,  V(f)  =  l/\R(f)\,  and  W(f)  =  1  (or  equivalently, 
U(f)  =  V(f)  =  | R(j)\V'  and  W(f)  =  1),  the  generalized 
correlation  spectrum  becomes: 

C(f  S(f)Rff)  (10) 

I  Rif)  I 

In  this  operation,  the  incoming  signal  is  correlated  with  a 
phase-only  version  of  the  replica,  hence  the  name  “phase- 
only.”  Since  the  incoming  signal  may  differ  from  the  replica 
greatly  in  magnitude,  the  POMF  may  perform  poorly  in  some 
cases. 

Eq.  (7)  can  be  viewed  as  a  combination  of  the  impulse 
response  of  Eq.  (9)  and  the  phase  only  matched  filter  of  Eq. 
(10)  as  a  function  of  the  replica  amplitude. 


(d)  Symmetric  Phase-Only  Matched  Filter  (SPOMF).  In 
the  final  case  where  U(f)  =  l/\S(f)\,  V(f)  =  l/\R(f)\,  and  W(f)  = 
7,  the  generalized  correlation  spectrum  becomes: 

S(f)Rff)  (11) 

I  S(f)  II  R(f)  I 

A  variation  of  this  filter  is  to  normalize  the  incoming  signal 
and  the  replica  spectra  with  the  square  root  of  their  respective 
magnitude  spectrum,  that  is,  U(f)  =  \S(f)[y\  V(f)  =  \R(f)\'v\ 
and  W(f)  =  1. 

In  Eq.  (11),  S(f)/\S(f)\  removes  the  magnitude  information 
from  the  incoming  signal  spectrum  and  retains  only  the  phase 
information.  Similarly,  R(f)/\R(f)\  removes  the  magnitude 
information  from  the  code  replica  spectrum  and  retains  only 
the  phase  information.  The  equalization  is  applied  to  both  the 
incoming  signal  and  replica,  hence  the  name  “symmetric.” 
Since  both  the  incoming  signal  and  replica  amplitude 
contents  are  involved  in  weighting,  the  spectral  filtering  is 
therefore  “balanced”  [Wernet,  2005]. 

Ideally,  a  flat  infinite  spectrum  produces  a  Dirac  delta 
function  in  the  time  domain.  By  eliminating  the  shape 
information  from  the  two  input  spectra,  the  phase  only 
filtering  (attenuating  the  magnitude  information  and 
accentuating  the  phase  information  in  the  frequency  domain) 
can  sharpen  the  correlation  peak  in  the  time  domain. 

The  importance  of  phase  in  signals  has  been  recognized  for 
many  applications  in  which  a  signal  can  be  recovered 
completely  or  in  part  from  knowledge  of  its  phase  alone.  In  a 
number  of  contexts  [9],  the  Fourier  transform’s  phase  data 
contain  more  of  the  “important”  information  than  the  Fourier 
transform’s  magnitude  data.  In  a  sense,  the  phase  reflects  the 
location  of  “events”  more  than  magnitude  whereas  the 
magnitude  contains  information  more  relevant  to  the  size  and 
shape  of  an  object.  The  time  shift  property  is  an  example:  a 
translation  in  position  (time  or  space)  of  a  signal  has  no  effect 
on  the  Fourier  transform  magnitude  but  only  affects  the  phase 
by  adding  a  linear  phase  term. 

Conditions  are  given  in  [9]  under  which  the  spectral 
magnitude  is  uniquely  specified  to  within  a  scaling  factor  by 
its  phase  function.  For  example,  the  log  magnitude  of  the 
Fourier  transform  is  the  Hilbert  transform  of  the  phase  of  a 
signal  with  all  poles  and  zeros  lying  only  in  the  left  half  or 
only  in  the  right  half  of  the  s-plane  (the  minimum-  or 
maximum-phase  condition). 

Since  the  autocorrelation  function  of  phase-only  signals  is 
always  an  impulse,  this  feature  has  been  used  in  designing 
methods  for  image  registration  and  recognition  [2]  and  digital 
image  velocimetry  [13]. 

IV.  Simulation  Results  and  Analysis 

Three  simulation  examples  are  provided  to  illustrate  the 
use  of  SPOMF  for  processing  GPS  signals  (BPSK  and  BOC- 
codes)  and  in  particular  its  superior  multipath  performance. 

A.  Bandwidth  Limiting  vs.  Phase-Only 

It  is  well  known  that  both  the  GPS  signal  and  the  RF  front- 
end  of  a  GPS  receiver  have  limited  bandwidth,  too.  The  effect 
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of  such  a  bandwidth  limitation  is  a  roimding  of  the  correlation 
function  at  the  peak.  To  see  this,  the  following  simulation  is 
conducted  in  which  the  C/A-code  of  GPS  SVN10  is  chosen 
and  the  sampling  rate  is  set  at  f  =  5  MHz  with  an  equivalent 
correlator  spacing  of  0.2046  chips.  A  5th-order  Butterworth 
low-pass  filter  is  used  for  band-limiting  filtering  with  its 
corner  frequency  set  at  1.023  MHz. 

Fig.  2  shows  the  first  150  samples  of  the  ideal  C/A-code 
(blue)  in  comparison  to  the  bandpass  filtered  code  (green). 
For  the  filtered  codes,  there  are  ringings  after  each  transition 
and  the  filtered  code  sequence  seems  to  be  shifted  in  time  due 
to  the  group  delay  impacted  by  the  lowpass  filtering.  The 
ringing  goes  away  and  transitions  of  first  order  appear  if  the 
corner  frequency  is  set  smaller  (or  the  transition  band  is  set 
larger).  Fig.  3  shows  the  ideal  code  amplitude  spectrum 
(blue)  and  the  filtered  code  spectrum  (green).  The  latter  is 
effectively  band-limited  around  1.023  MHz  as  designed.  Fig. 
4  shows  the  ideal  correlation  function  (blue)  of  a  triangular 
shape  and  the  filtered  correlation  function  (green),  which  is 
smaller  in  amplitude  and  is  rounded  up  at  the  correlation 
peak.  Fig.  5  compares  the  autocorrelation  functions  of  the 
ideal  PRN  code  (blue),  the  filtered  code  (green),  the  phase- 
only  version  of  the  ideal  code  (red),  and  the  phase-only 
version  of  the  filtered  code  (cyan).  The  autocorrelation 
functions  of  the  two  phase-only  codes  have  the  same  shape, 
which  is  much  narrower  and  of  smaller  noise  floor  than  the 
two  originals. 

To  visualize  a  PRN  code  and  its  phase-only  version,  the 
following  simulations  are  conducted  first  for  a  BPSK  code 
and  then  for  a  BOC  code.  Fig.  6  shows  the  GPS  C/A-code 
sampled  at  5  MHz  and  its  phase-only  version  (normalized  to 
unity  amplitude).  Their  normalized  autocorrelation  functions 
are  shown  in  Fig.  7.  Similarly,  Fig.  8  shows  the  GPS  C/A- 
code  sampled  at  10  MHz  and  its  phase-only  version 
(normalized)  with  their  normalized  autocorrelation  functions 
in  Fig.  9. 

Fig.  10  shows  a  BOC(10,  5)-code  sampled  at  50  MHz  and 
its  phase-only  version  (normalized  to  the  unity  amplitude). 
Their  normalized  autocorrelation  functions  are  shown  in  Fig. 
11.  Similarly,  Fig.  12  shows  the  same  BOC  code  sampled  at 
100  MFIz  and  its  normalized  phase-only  version  with  their 
normalized  autocorrelation  functions  in  Fig.  13. 

It  is  interesting  to  see  that  the  phase-only  version  of  both 
the  BPSK  and  BOC  codes  produces  a  pair  of  spikes  of 
opposite  directions  at  a  transition.  The  signal  returns  to  zero 
during  quiescent  periods  with  small  variations.  The  effect  is 
more  pronounced  when  the  sampling  rate  is  higher.  This 
indicates  that  the  phase-only  waveform  possesses  the 
property  of  some  specially  designed  codes  such  as  double 
delta,  strobe,  pulsed  aperture,  and  gated  correlator  that  are 
employed  in  the  state  of  the  art  conventional  GPS  receivers 
and  are  known  to  reduce  multipath  errors  [21,  22,  23]. 

B.  BOC  Modulation  Codes 

In  this  simulation,  we  use  the  following  method  to 
construct  a  BOC(10,  5)-code  similar  to  the  M-code.  It  is  easy 
to  generate  the  subcarrier  with  a  square  wave  at  the 
fundamental  frequency  of  10.23  MHz  (an  equivalent  chipping 


rate  is  20.46  Mcps).  We  use  the  CL  component  of  the  L2C 
code  (which  is  1.5  s  long)  as  the  PRN  code  at  5.115  Mcps 
(another  choice  is  the  P-code).  The  correlation  and  spectral 
properties  of  these  surrogate  codes  may  not  be  as  good  as  the 
actual  M-code,  which  is  infinitely  long  and  cryptographic, 
thus  representing  a  worse  case  of  signal  conditions. 

The  PRN-code  and  the  square  wave  are  modulo-2  added, 
multiplied  with  a  complex  exponential  to  simulate  the 
residual  carrier  plus  complex  Gaussian  white  noise.  The 
resulting  signal  at  the  baseband  is  then  sampled  at  50  MHz 
and  passed  through  a  5th-order  Butterworth  low-pass  filter 
with  its  corner  frequency  set  at  12  MHz.  Fig.  14  shows  the 
amplitude  spectrum  of  the  code  sequence  sampled  at  50 
MHz.  Fig.  15  shows  the  spectrum  after  the  signal  is  lowpass- 
filtered  within  12  MHz. 

In  this  simulation,  the  residual  Doppler  error  is  set  to  be 
200  Hz.  The  code  phase  error  is  14  the  sampling  interval  (i.e., 
between  two  samples).  The  initial  carrier  phase  is  drawn 
uniformly  from  [0,  27t).  A  complex  Gaussian  noise  of  unity 
variance  is  added.  The  signal  amplitude  is  adjusted  to 
simulate  the  desired  SNR  level  according  to: 

Cl  No  | 

a  =  (t,io  ■»  y  (12) 

where  7}  =  0.001  s  and  C/No  =  30  dB-Hz  for  this  simulation. 
The  resulting  SNR  is  1  (in  ratio). 

A  segment  of  1  ms  code  is  taken  from  a  long  sequence  of 
the  simulated  M-code  starting  from  the  1000th  chip  as  the 
replica  (the  true  peak  location  is  at  2444.3).  Fig.  16  shows  a 
portion  of  the  correlation  where  the  peak  value  is  about 
31200  and  the  noise  floor  peak  is  about  1900.  The  peak  to 
noise  peak  ratio  is  16.42  (or  24.31  dB). 

The  ideal  correlation  value  is  50000.  However,  the  code 
phase  error  of  14  the  sampling  interval  (0.2046  chips)  reduces 
it  by  a  factor  of  0.7954.  The  Doppler  error  of  200  Hz  (relative 
to  the  1  ms  integration  interval)  introduces  a  loss  factor  less 
than  0.92.  The  12  MHz  lowpass  filter  introduces  a  loss  factor 
less  than  0.86.  The  practical  peak  value  is  around  31466, 
which  is  close  to  what  we  can  observe  from  Figs.  16  and  17 
(the  absolute  peak  value  around  31200). 

It  is  also  clear  from  Figs.  16  and  17  that  the  BOC 
modulation  has  sidelobes  in  addition  to  the  main  peak.  This  is 
problematic  for  conventional  tracking  loops  with  the  risk  of 
locking  onto  a  secondary  peak  if  the  correlator  spacing  is 
small,  less  than  a  !4  of  the  subcarrier  wavelength  for  instance. 

Fig.  17  shows  the  real  and  imaginary  components  of  the 
complex  correlation.  In  the  acquisition  mode,  carrier  phase 
and  frequency  tracking  are  not  yet  applied.  It  is  not 
dominated  by  the  real  component  as  is  caused  by  the  Doppler 
frequency  error.  For  this  particular  segment  of  data,  the  real 
and  imaginary  components  have  about  the  same  power. 

Figs.  18  and  19  compare  the  conventional  (green) 
correlation  with  the  impulse  response  (red)  and  the  symmetric 
phase-only  (blue)  correlation.  As  explained  before,  similar  to 
the  phase-only  correlation  in  Eq.  (10),  the  impulse  response 
in  Eq.  (9)  equalizes  the  replica  spectrum.  Although  producing 
very  narrow  peak,  it  tends  to  develop  large  sidelobes  as 
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Fig.  2.  Ideal  PRN  vs.  Lowpass  Filtered  Code 


5th-butterworth  lowpass  filter,  con=1.023MHz,  fs=5MHz 


Fig.  4.  Ideal  PRN  vs.  Lowpass  Filtered  Autocorrelations 
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Fig.6.Ideal  PRN  vs.  Phase-Only  Codes 
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Fig.  3.  Ideal  PRN  vs.  Lowpass  Filtered  Spectra 


5th-butterworth  lowpass  filter,  con=1.023MHz,  fg=5MHz 


Fig.  5.  Conventional  vs.  Phase-Only  Autocorrelations 


normalized  peak  for  BPSK  code  vs.  phase-only  code 


Fig. 7. Ideal  vs.  Phase-Only  Autocorrelations 
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code  sequenece  vs.  normalized  phase-only  code  normalized  peak  for  BPSK  code  vs.  phase-only  code 


Fig.8.  Ideal  PRN  vs.  Phase-Only  Codes  Fig-9.  Ideal  vs.  Phase-Only  Autocorrelations 
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Fig.  10.  Ideal  PRN  vs.  Phase-Only  Codes 


Fig.  1 1 .  Ideal  vs.  Phase-Only  Autocorrelations 
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Fig.  12.  Ideal  BOC(10,  5)  vs.  Phase-Only  Codes 


Fig.  13.  Ideal  vs.  Phase-Only  Autocorrelations 
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shown  in  Fig.  18.  Spectral  filters  such  as  zone  zeroing  and 
individual  excision  have  been  used  to  reduce  the  noise  floor 
[18]. 

The  symmetric  phase-only  matched  filter  equalizes  both 
the  incoming  signal  and  replica  spectra.  As  such,  it  not  only 
produces  a  peak  as  narrow  as  the  phase-only  matched  filter, 
blit  also  keeps  the  noise  floor  as  low  as  the  conventional 
correlator  as  shown  in  Fig.  19. 

C.  BPSK  Modulation  Codes 

The  multipath  performance  of  the  SPOMF  is  demonstrated 
in  the  following  example  with  a  BPSK  code  (the  GPS  SVN10 
C/A-code).  The  signal  ms  boundary  is  ahead  of  the  first 
sample  by  0.4 5TS.  One  multipath  component  is  considered 
with  T/  =  0.5  chips  (2.4438  samples)  and  =  0.2. 

We  first  compare  the  signal  channel  impulse  response  with 
the  conventional  correlation  function  in  Figs.  20  and  21  for 
two  spectral  filtering  methods,  namely,  zone  zeroing  and 
individual  excision  [17,  18].  The  top  plots  in  Figs.  20  and  21 
correspond  to  the  case  where  the  incoming  signal  contains 
direct  signal  and  noise  only  without  multipath  whereas  the 
bottom  plots  of  Figs.  20  and  21  contain  the  direct  signal, 
noise,  and  multipath. 

It  is  clear  that  the  peak  of  the  impulse  response  function  is 
narrower  than  that  of  the  correlation  function.  The  latter  is 
about  2x5/1.023  =  9.7752  samples  wide  whereas  the  former 
is  about  3  samples  wide.  This  has  significant  ramifications  in 
multipath  performance. 

For  multipath  delays  larger  than  one  sampling  interval, 
which  is  0.2  chips  in  this  simulation,  the  individual  multipath 
components  can  be  resolved,  thus  not  causing  errors  to  the 
direct  signal.  For  those  within  0.2  chips,  the  multipath  affects 
the  correlation  values  at  adjacent  sampling  points.  Since  the 
direct  signal  is  stronger,  an  interpolation  can  be  used  to  find 
the  underlying  peak  location. 

In  Fig.  22,  the  multipath  delay  is  varied  from  0.05  to  1.5 
chips  in  step  of  0.1  chips.  In  each  case  we  show  the 
performance  of  four  different  methods.  The  four  methods  are: 
(1)  quadratic  fitting  to  the  correlation  function  (blue),  (2)  sine 
fitting  to  the  impulse  response  without  filtering  (green),  (3) 
sine  fitting  to  the  impulse  response  with  zone  zeroing  (red), 
and  (4)  sine  fitting  to  the  impulse  response  with  individual 
excision  (cyan). 

The  correlation  function  interpolated  estimate  without 
multipath  mitigation  follows  closely  the  multipath  error 
envelope,  which  is  determined  by  the  correlation  spacing  d  = 
0.2  chips  and  multipath  strength  Oy  =  0.2  (i.e.,  da )  =  0.04  or 
12  m). 

The  two  spectral  filtering  methods  perform  consistently 
better  than  the  one  without  spectral  filtering.  Between  the  two 
spectral  filtering  methods,  the  individual  excision  method 
seems  marginally  better  than  the  zone  zeroing.  The  data 
shown  in  the  figures  are  taken  from  sample  runs,  each  with  1 
ms  worth  of  data  samples.  Averaging  over  many  runs  may 
smooth  out  variations  in  the  curves. 

In  Fig.  23,  the  delay  estimation  errors  of  the  symmetric 
phase-only  matched  filter  are  compared  with  those  of  two 
spectral  filtered  impulse  responses.  Once  again,  the 


individual  excision  method  marginally  performs  better  than 
the  zone  zeroing.  The  symmetric  phase-only  matched  filter 
performs  best. 

In  addition  to  the  sample  behavior  shown  above,  we  also 
conducted  Monte  Carlo  runs.  The  root  mean  squared  (RMS) 
values  of  the  errors  in  estimating  the  direct  signal  delay  (z0) 
as  a  function  of  the  multipath  signal  delay  (zf)  are  shown  in 
Figs.  24  and  25  for  C/N0  =  40  dB-Hz  and  30  dB-Hz, 
respectively.  The  multipath  strength  is  fixed  at  at  =  0.2. 

In  each  figure,  there  are  three  pairs  of  curves,  representing 
three  different  estimation  techniques.  Each  pair  consists  of 
two  types  of  signals,  namely,  (1)  direct  signal  plus  noise  (s  + 
n)  and  (2)  direct  signal  plus  noise  plus  multipath  signal  (s  +  n 
+  m).  The  first  delay  estimation  technique  is  the  normalized 
early  minus  late  delay  error  discriminator  applied  to  the 
correlation  power.  With  the  correlation  spacing  d  =  0.2046 
chips  and  the  multipath  strength  at  =  0.2 ,  the  expected  error 
is  on  the  order  of  daj  =  0.04  or  12.28  m. 

The  second  technique  is  the  fitting  of  a  sine-function  to  the 
impulse  response  peak  after  individual  excision  of  spectral 
spikes  is  done  in  the  frequency  domain.  The  third  technique 
is  the  fitting  of  a  sinc-frinction  to  the  symmetric  phase-only 
matching  peak. 

Each  Monte  Carlo  simulation  consists  of  100  runs.  For 
each  run,  the  initial  carrier  phase  and  noise  are  drawn 
randomly  from  their  respective  distributions.  However,  the 
random  values  are  kept  the  same  when  the  multipath  signal 
delay  is  varied  from  0.1  to  1.4  chips  in  the  step  of  0.1  chips. 
This  helps  explaining  why  the  delay  estimation  errors  for  the 
three  “signal  plus  noise”  cases  (green,  cyan,  and  yellow, 
respectively)  remain  constant  in  the  figures. 

While  a  conventional  correlation-based  delay  estimation 
technique  produces  the  typical  multipath  error  envelope  even 
when  the  SNR  is  high  (see  the  blue  curve  of  Fig.  24),  both  the 
impulse  response  (red)  and  symmetric  phase-only  matching 
(purple)  techniques  perform  similarly  with  an  almost  flat 
error  level,  indicating  the  insensitivity  to  multipath.  Their 
error  level  is  slightly  higher  than  the  “signal  plus  noise” 
cases.  As  expected,  they  all  attain  the  same  level  of 
performance  when  the  multipath  delay  goes  beyond  (1+d) 
chips  after  which  the  multipath  does  not  interfere  with  the 
direct  signal  and  the  noise  dominates. 

When  the  SNR  is  low,  all  “signal  plus  noise”  curves  are 
raised  in  their  error  level.  The  impulse  response  (red)  and 
symmetric  phase-only  matching  (purple)  techniques  start  to 
behave  differently.  Both  the  techniques  are  still  effective  in 
suppressing  the  multipath  effects  because  their  “signal  plus 
noise  plus  multipath”  curves  (red  and  purple)  stick  rather 
closely  to  the  “signal  plus  noise”  curves  (cyan  and  yellow). 
However,  the  impulse  response  (red  and  cyan)  technique 
experiences  a  noise  amplification  effect  as  observed  in  the 
sample  behavior.  The  error  level  (red  and  cyan)  is  higher  than 
the  conventional  correlation  noise  floor  (green). 

In  contrast,  the  symmetric  phase-only  matching  (purple 
and  yellow)  technique  maintains  an  almost  flat  multipath 
error  level,  which  is  even  slightly  lower  than  the  conventional 
correlation  noise  floor  (green).  This  may  be  explained  by  the 
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"m-code"  sampled  at  50  MHz,  full  band 


Fig.  14.  M-Code  Sampled  at  50  MHz 
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Fig.  15.  Spectrum  after  Filtering 
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Fig.  16.  A  Portion  of  Correlation  at  50  MHz 
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Fig.  17.  Real  and  Imaginary  Parts  of  Correlation  (50  MHz) 
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generalized  frequency-domain  correlator  (GFDC),  C/NQ=30dB-Hz 


Fig.  18.  Conventional  vs.  Phase-Only  Correlations 


Fig.  19.  Conventional  vs.  Phase-Only  Correlations 


195 

9 


delay  estimation  rms  error,  m  timing  error, 


signal  &  noise 


signal,  multipath  &  noise 


delay,  0.2046  chips 


Fig.  20.  Ideal  vs.  Phase-Only  Correlation 
(Impulse  Response  with  Zone  Zeroing) 
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Fig.  21.  Ideal  vs.  Phase-Only  Correlation 
(Impulse  Response  with  Individual  Excision) 


multipath  insensitive:  ifft[FS./FR] 


GFDC:  multipath  desensitive  delay  estimation 


"ig.  22.  Multipath  Error  Envelopes  for  Ideal  vs.  Phase-Only 
Correlation  with  Different  Filtering 
(Impulse  Response  0Ci  =  0.2,  C/No  =  30  dB-Hz) 


Fig.  23.  Multipath  Error  Envelope  for  Filtered  Phase-Only  vs. 
Symmetric  Phase  Only  Correlation 
(oci  =  0.2,  C/N0  =  30  dB-Hz) 


monte  carlo  simulation:  a=0.2,  C/N0=40dB-Hz  ,100  runs  monte  carlo  simulation:  a=0.2,  C/NQ=30dB-Hz  ,100  runs 


Fig.  24.  Delay  Estimation  Error  RMS 
(a[  =  0.2,  C/No  =  40  dB-Hz) 


Fig.  25.  Delay  Estimation  Error  RMS 
(a[  =  0.2,  C/N0  =  30  dB-Hz) 
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fact  that  the  slope  of  SPOMF  peak  is  steeper  than  that  of  the 
conventional  correlators,  thus  having  smaller  noise  variance. 
Overall,  the  simulation  results  indicated  the  consistency  and 
superiority  of  the  SPOMF  technology  versus  conventional 
correlators. 

V.  Conclusions 

In  this  paper,  we  set  forth  a  generalized  frequency-domain 
correlator  (GFDC)  as  an  efficient  computing  engine  for  a  new 
generation  of  software  GNSS  receivers.  With  its  inherent 
flexibility,  the  GFDC  can  be  programmed  to  provide  the  most 
suitable  form  of  correlation  when  operating  in  both  the 
acquisition  and  tracking  modes.  In  acquisition,  the  desired 
correlation  is  of  a  wide  base,  which  can  be  obtained  by  FFT- 
implemented  conventional  correlation  for  the  BPSK 
modulation  or  for  a  single  sideband  if  the  BOC  modulation. 
In  tracking,  the  desired  correlation  is  of  a  narrow  base,  which 
is  achieved  with  the  symmetric  phase-only  matched  filtering 
(SPOMF).  The  same  operation  is  applicable  to  both  the 
BPSK  and  BOC  codes.  With  spectral  filtering  and  feedback 
closure,  the  GFDC  realizes  a  seamless  transition  from 
acquisition  and  tracking. 

Simulation  results  showed  the  SPOMF  indeed  provided  a 
very  sharp  correlation  peak  that  was  more  accurate  in  timing 
and  less  sensitive  to  multipath.  The  SPOMF  has  been  applied 
to  real  GPS  data  that  are  known  to  contain  multipath.  The 
initial  test  results  will  be  presented  in  [20]  together  with 
analysis  of  such  effects  as  signal  bandwidth  and  code  replica 
filtering. 
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